PEPS as unique ground states of local Hamiltonians 
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In this paper we consider projected entangled pair states (PEPS) on arbitrary lattices. We 
construct local parent Hamiltonians for each PEPS and isolate a condition under which the state is 
the unique ground state of the Hamiltonian. This condition, verified by generic PEPS and examples 
like the AKLT model, is an injective relation between the boundary and the bulk of any local region. 
While it implies the existence of an energy gap in the ID case we will show that in certain cases 
(e.g., on a 2D hexagonal lattice) the parent Hamiltonian can be gapless with a critical ground state. 
To show this we invoke a mapping between classical and quantum models and prove that in these 
cases the injectivity relation between boundary and bulk solely depends on the lattice geometry. 



I. INTRODUCTION 

Exactly solvable models for quantum mechanical spin 
systems are important testbeds and sources of insight for 
the theory of condensed matter systems and their statis- 
tics. Unfortunately, such models are rare even in ID 
and, disregarding free particles, essentially non-existent 
in higher spatial dimensions. If the focus lies in low tem- 
perature physics then quasi-exactly solvable models for 
which at least the ground state is known exactly become 
interesting. A powerful tool for constructing such mod- 
els in ID is the matrix product state (MPS) formalism 
P, [2, [1| which has his roots in the AKLT model [1] and 
has experienced a remarkable development in the last 
years concerning both analytical results Q and numer- 
ical methods @. In two and higher spatial dimensions 
a generalization of the MPS idea is known as projected 
entangled pair state (PEPS) formalism 0, H, 0], earlier 
attempts of which can already be found in the seminal 
work of AKLT 3 as well as in studies of dimer models 
[lo| and antiferromagnetic vertex spin models [11,] . The 
idea of the PEPS representation of quantum states has 
its roots in quantum information theory but it found 
most of its app lications in quantum many-body physics 

0J,i[ll,[li. 

The present work will continue on this line and investi- 
gate quasi-exactly solvable models based on PEPS. Our 
main interest lies hereby in the uniqueness of the ground 
state. We prove that 

• every PEPS satisfying a condition we call 'injectiv- 
ity' is the unique ground state of a local Hamilto- 
nian. 

• For quantum states which are coherent versions of 
classical Gibbs states we show that injectivity de- 
pends solely on the lattice geometry. 

• By invoking the mentioned classical-to-quantum 
mapping [8[ we construct a PEPS which is the 



unique critical ground state of a local Hamiltonian 
on a 2D hexagonal lattice. This shows that, in con- 
trast to ID MPS, injectivity does not imply the 
existence of an energy gap. 

• We give an example on the 2D square lattice where 
uniqueness can be proven in the absence of injec- 
tivity. 

• We provide a computable sufficient condition for 
the existence of an energy gap. 

• In the appendix we show that every translational 
invariant PEPS on a 2D square lattice admits a 
representation with site-independent tensors. 

Note that the relation between injectivity and unique- 
ness generalizes in a fairly direct way what is known in 
the ID case [l|, to arbitrary lattice geometries. The 
main difference lies in the fact that beyond the ID case 
the existence of a gap above the ground state energy be- 
comes a rather independent property, the determination 
of which essentially remains an open problem. 



II. PRELIMINARIES 

Let us start by recalling that a PEPS is a quantum 
state \(p) constructed over a graph (Fig.l) such that 
each vertex corresponds to a physical site and the ge- 
ometry of the graph typically resembles the spatial ge- 
ometry of a lattice. To construct the PEPS we assign 
to each edge (bond) of the graph a maximally entan- 
gled state J2iLi N)- this way, each vertex gets as- 
signed to as many D-dimensional systems as there are 
adjacent edges — we may think of these as a virtual sub- 
structure. To obtain the final state we apply a map 
P : (8) • • • (g) — > at each vertex which maps 
the virtual onto the physical system. Note that P and 
d may depend on the vertex and D on the edge. For 
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the sake of simplicity we will, however, choose the di- 
mensions D and d to be constant over the entire graph. 
Let us denote by the number of edges at vertex v and 
parameterize the corresponding map P^^^ by 

pW=^ E A(5|z)(ji,...,JeJ 

such that each A^^^ is a tensor with indices j — 
(ji, . . . , je^). From this we obtain a PEPS of the form 



(1) 



11,12, ... = 1 



where C means the contraction of all tensors A'f^ accord- 
ing to the edges of the graph. If, for instance, we deal 
with an iV X M square lattice on a torus we may specify 
each vertex by the respective row j and column k such 
that 



i^> - E c 



(2) 

where each A^^'"^' is now a 4-index tensor which is con- 
tracted according to the lattice: 



l(l,2) 

*(1,2) 



A 



(2,1) 

*H2,1) 



(2,2) 

*i(2,2) 



Clearly, if the tensors A)^' do not depend on the ver- 
tex V = (j, k) then the corresponding PEPS \ f) is trans- 
lational invariant. Conversely, if \lp) is translational in- 
variant, then there exists always a representation (with 

presumably larger D) such that A'f'^^ = Ai is site- 
independent. A proof of this statement can be found 
in the appendix. We will throughout not assume such a 
symmetry and allow the Aj's to be site dependent though 
we typically omit to write the explicit dependence A. 



(v) 





FIG. 1: Left: A PEPS can be assigned to an arbitrary graph. 
Each vertex corresponds to a physical site and every edge 
depicts the contraction of tensors assigned to the vertices. 
Right: A parent Hamiltonian which has the PEPS as unique 
ground state can be constructed by first joining vertices in 
injective regions (now surrounded by solid lines) and defining 
interactions (depicted by bars) between connected neighbor- 
ing regions. 



which allows us to construct parent Hamiltonians for ev- 
ery PEPS \ip). Consider a region R on the lattice with re- 
spective reduced density operator = iY\^\Lp){Lp\. The 
support space Sr of pr has dimension at most Z?*^" where 
cr is the number of edges connecting the interior of R to 
the exterior. Denote by \R\ the number of lattice sites in 
R with total Hilbert space TLr = 'C<^®W . Since Sr grows 
as the boundary of R whereas T-Lr grows like its volume, 
we can always find a sufficiently large region such that 
PR will have a non-trivial kernel ker(pfl') — T-Lr \ Sr. It 
is thus a ground state of every Hamiltonian h with the 
following properties 

• ft. > 0, 

• ker(/i) — Sr, 

since by construction tr [pR.h] = 0. This Hamiltonian can 
then be extended to the whole lattice by assigning such 
a surrounding region R^ and a corresponding to every 
vertex v and defining 



H 



E^- 



(3) 



In this way we have again H > and {(p\H\if) = 
J2v^^iPRv^v] = 0. The Hamiltonian is thus local and 
frustration free in the sense that \^p) minimizes the 

energy locally for every interaction term . The state is 
therefore sometimes called an optimum ground state [ll[ . 
We remark that it was recently proven in [14] that all lo- 
cal gapped Hamiltonians can be efficiently approximated 
by frustration free ones. 



III. PARENT HAMILTONIANS 



IV. INJECTIVITY 



In the following we will mainly consider lattices instead 
of arbitrary graphs as this gives us a meaningful notion 
of locality and ensures that volumes grow faster than 
boundary areas. The latter is, in fact, the main property 



We will now discuss the main condition which will later 
on turn out to be sufficient for proving uniqueness of the 
ground state for Hamiltonians of the form in Eq. ^ . Con- 
sider a PEPS \(f) given by tensors Ai (depending on the 
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vertex). Let i? be a connected region on the graph con- 
taining \R\ sites (vertices) and having outgoing edges 

(bonds). The hnear map Tr ■■ C^'" ^ defined by 



(4) 



maps the virtual boundary of R to the physical bulk. 
Here stands for the tensor obtained by contracting all 
the Ai 's in the region R with bonds open in the boundary 
of R (and then C is a tensor with the same number en 
of bonds). Graphically we will represent C[v4ij,C] as 



C 



where we have joined all the cr open bonds of Ai^ and 
C in one line. 

If Tr is an injective map, we will simply say that the 
region R is injective. Similarly, if there is a proper cov- 
ering of disjoint injective regions for the entire graph we 
will call the PEPS itself injective. 

Of course, the injectivity property of R is equivalent to 
the fact that the finite set {Ai^}i^ generates the space 
of tensors with e/j ZJ-dimensional indices. Just counting 
degrees of freedom one can see that the dimension of the 
space of tensors that we have to generate grows like the 
boundary D'^'^, whereas the number of tensors that we 
have grows like the volume d'^L This strongly suggest 
that the injectivity condition is generic on any lattice, 
i.e., almost always fulfilled as long as d'^' > D'^'^ . Indeed, 
if one takes a numerically random PEPS, it appears to be 
fulfilled in all the cases. As a more concrete example, one 
can also check that the 2D AKLT state [4J| is injective. 

In the remaining part of this section we will prove a 
pair of Lemmas which are particularly concerned with 
lifting the injectivity property and its consequences from 
smaller to larger regions. The reader not interested in 
the technical details may skip this part. 

Lemma 1. If R and S are two disjoint injective regions, 
then RU S is also injective. 

Proof. We start from the picture 



A, 



A, 



Since both Ai^^ and Ai^ generate the respective space 
of tensors, then, summing in i/j and is we can get 



for each ao and cq. Hence, we have injectivity in the 
region R U S, since we can generate all the canonical 
basis vectors of the respective space of tensors. □ 



Let us call Gr = range(ri{) the subspace generated 
in the physical bulk by acting on the virtual boundary. 
Note that the reduced density operator pR is supported 
in Gr 3 Sr. Moreover, if R and its complement are 
injective then the converse holds as well so that Gr — Sr. 

Lemma 2. Let Ri, R2, R3 be three disjoint regions, and 
Ti-R — C*' ' denote the total Hilbert space corresponding 
to region R. 

1- Gr^uR2UR3 ^ {GR1UR2 ® Ti-Rs) n {TCri (8) Gr^uRs)- 

2. If Ri and R3 are not connected and R2 
and i?3 are both injective, then Gr-^uR2UR3 = 
(G_RiUfl2 T^Rs) 1^ C^Ri Gr^uRs)- 

3. If all three regions are injective, then 

GR1UR2UR3 = {GriUR2 ^ T^R3)^{'^Ri ® Gr.^ur^)P\ 
{GriURs ® T^Ra)- 

Proof. It is clear that € Gr^ijr^ ® T^Rs if and 
only if lip) can be written as J2iii^^RiuR2{GiR^)\iR3)- 
Then, if 1-0) G Gr-^[jR2UR3, since we can write \ip) — 
T.tn,,tR^,iR^ C{AR^AR^AR^X)\iR^iR^iR^), calling C,^^ = 
C{Ai^^X), we have that £ Gr,uR2 ® T^Rs- 

Of course, the same can be done for Hr^ (X" GR^uRa, 
which proves part 1. of the Lemma. 

For the other parts note that if IV') <= 
{Gr,uR2 T^Rs) ^ C^Ri ® Gr^ur^), then 

IV^) ^^^RlUR2{D,^J\iR^) = ^rfl,,ufl,3(C»Bi)l*fl.l). 

which imphes that 



^iR2 


d 




A, 


R3 



A, 




U 




A, 


«2 



for every iR^ , iR.^ , iR,^ . Or equivalently (summing in the 
repeated indices), 

j^,l.d j^d,x,y ^u,l^x.y j^,h,u j^,l^d j-ja,h,l^d ^f^^ 



'I.R2 



where the sum over a and y is only non-trivial if i?i is 
not separated from R^. 

In order to prove 2. we now use injectivity in R2 and 
i?3. That is, for every fixed set of indices lo,UQ,dQ and 



and 



IR3 



such that 



xq, yo there exist a, 

E uo,lo,do AU,l,d _ s- XX 
^in^ - du,uoOl,loOd,do, 



xo,yo Ad,x,y 



(6) 
(7) 
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Inserting this into Eq. (fT5|) we get that 



(note that M"-''''''^'^ does not depend on d). Hence, the 
coefBcients of '0 have the form of the left part of the 
foUowing identity (proven below) 





(8) 



which implies that under the conditions of part 2. S 
G_RiUiJ2Ui?3 completing the proof of the second part. 

For 3. we first repeat the above proof with a cyclic per- 
mutation of the regions, i.e., {Ri, R2, R3) (i?2, ^3, Ri) 
leading to the identity in ([5]). Then we use injec- 
tivity in all three regions in order to replace A^^'^, 

Ktf ^'^^ Kr.^ ^a^aah.ha-. ^i-h ^^d (5j;,2;o'5iy,yo respec- 
tively. Inserting this into Eq.® we obtain M"-''''^'^'^ — 
Ja.v^f"''''^'^'" so that 



where the sum runs over all edges (a, /3) in the super- 
lattice and ha, 13 > is any nearest- neighbor interaction 
term satisfying keT{ha,[3) = GR^uR,i- 

Theorem 3. For every injective PEPS on an arbitrary 
lattice there is a local frustration-free Hamiltonian such 
that the state is its unique ground state. 

Remark: Note that though Eq.Q gives a standard 
construction this Hamiltonian is not unique. We may in 
particular choose h'^ ^ = ha,f}-\-ha,f}Qha,p for any Q > 0. 
Similarly, we can use the same construction for larger 
regions, or sometimes even for smaller ones for which 
injectivity does not yet hold (see below). The advantage 
of the latter is to reduce the number of sites involved in 
the interaction, as done for instance in the AKLT model. 

Proof. Using the Hamiltonian of Eq.® and writing, with 
a slight abuse of notation, Gr for Gr ® 'H\b we have to 
show that 



G, 



(10) 



as the left hand side is the ground state space of H and 
the right hand side is the one-dimensional space corre- 
sponding to the PEPS. We will show this identity by 
induction. To this end consider a collection of disjoint 
regions LifLf^Ai — A such that Ai, . . . , A^ are connected 
to a region B (all regions being injective). Then as 
Ga ^ GaouA, we have that 



,«i?2 '^-^3 




\^Rl,^R2^^R3 



is indeed in Gijiufl2Ufl3- □ 



V. UNIQUENESS OF THE GROUND STATE 

Our aim now is to prove that every injective \(p) is 
the unique ground state of a Hamiltonian of the form 
in Eq.([3]). Consider an injective tile of the lattice, i.e., 
a covering by disjoint injective regions and let us merge 
sites within each injective region such that in the new, re- 
grouped lattice each single vertex is injective. We connect 
the vertices in this super-lattice if the respective regions 
were connected before (see Fig.l). As a Hamiltonian we 
choose 



H 



E 

(",/3) 



ha,p > 



(9) 



N N N 

Ga Pi GAiUB C P (^GAoUAit^GAiUB^ — P GaoUA.UB, 
i=X i—1 i—1 

(11) 

where the last equality follows from part 2. of Lemma 
[21 If we exploit in addition that Ga Q Gm^iAi for any 
/ C {0, . . . , N} together with part 3. of Lemma [2] we 
obtain 



N 



Ga P GAiUB c Gaub- 



(12) 



4=1 



Eq. (jl2p can now be used as an induction step in order to 
show that the l.h.s. of Eq. fTU)) is contained in the r.h.s. . 
The converse inclusion follows, however, from part 1. of 
Lemma [5] which completes the proof. □ 



VI. QUANTUM STATES FROM CLASSICAL 
MODELS 



In y we constructed PEPS from any classical spin 
model with nearest-neighbor interaction (with d possible 
configurations per site and a fixed given temperature ^) 
in such a way that (i) the physical dimension coincides 
with the bond dimension, i.e., d = D; and (ii) the PEPS 
reproduces the correlations of the classical thermal state. 
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In a sense, this map replaces thermal by quantum fluc- 
tuations. In this section we will analyze in which cases 
PEPS constructed in this way are unique ground states 
of their parent Hamiltonians from Eq. ^ . 

First of all we will review the construction of 
For a classical Hamiltonian of the form H{x) = 
'^(^^ jjh{xi,Xj), with Xi ~ we define the asso- 

ciate PEPS as 





H) = cxp 



and hij is a diagonal operator 
that acts as hij\xiXj) = h(xi,Xj)\xiXj). To see explicitly 



where |+) = '£x=i \^) 

the PEPS structure of {ip) one can notice that the non- 
unitary gate exp[— ■l^.y] is indeed equal to Pi ® Pj\I), 
where |/) — l-'^^) auxiliary maximally entan- 

gled state between the sites i and j and Pi : — > 
is an operator (from the joint physical-virtual system 
into the physical system) acting as Pi\xiv) = \xi){(f>l..\v). 
Here the vectors \4>l..) can be obtained by a SVD to verify 
eyip[—^h{xi,Xj)] = J^ki'l^xil^) i^l^xj) ■ Then, the tensors 
defining the PEPS \^) are 



Al 



(13) 



where |(/)|.) may, of course, depend on the lattice site 
z, and the product is taken over all edges e connected 
with this site. Generically, c/)'^ can be considered (via 
{xi\ I— > {4>%.\) an invertible d x d matrix. For instance, 
for the Ising model H{x) — — j) ^i^jj — =^1, 



sinh ^ 



— A/sinh ^ 



which is indeed invertible for all /3. 



cosh ^ 
cosh ^ 



FIG. 2: Examples of 2D lattice geometries leading to injectiv- 
ity when the PEPS is constructed from a classical spin model 
with corresponding interaction graph. A 3D example would 
be a tetrahedral crystal lattice (ice) . 



where the bar always means "at the boundary" and C{x) 
is the result of the contractions along all e ^ E. In the 
case \E\ = \R\ we have that the vectors x and a have the 
same dimension and F = •SieeEff''^- 

Hence, F is invertible so that Tji is injective [l^. On 
the other hand if |£'| > \R\ then F is rectangular and 
not invertible, such that Tji cannot be injective. Hence, 
we do not have injectivity for the square lattice but we 
do have it e.g. for the tetrahedral 3D or the hexagonal 
2D lattice (Figure Oa). A 2D square lattice can as well 
lead to injectivity if we either introduce a substructure at 
each point (Figure [He) or if we allow for defects (Figure 
[2]b) : if with probability p a bond is missing then injective 
regions contain of the order of (1 —pf'/p^ lattice sites. 



With this at hand we can now give an example of an 
injective PEPS whose associated parent Hamiltonian is 
gapless. It is the PEPS associated to the 2D classical 
(isotropic) hexagonal Ising model. On the one hand, we 
have injectivity and hence the PEPS is the unique ground 
state of its parent Hamiltonian. On the other hand, the 
Ising model becomes critical for /3 = i ln(2 -f \/3) [l5|, 
which implies that the PEPS has correlations with power- 
law decay. By [l^ , this implies that the parent Hamilto- 
nian has to be gapless. 



A. Injectivity and Criticality 



Surprisingly injectivity does, in the case of PEPS cor- 
responding to classical models, not depend on the type 
of interaction, but merely on the lattice geometry of the 
PEPS, which coincides with the interaction graph of the 
classical model. If a subset on the interaction graph is 
such that every site has at most one outgoing edge, then 
the corresponding region has the injectivity property — 
otherwise it does not. The simple reason is as follows. 

Consider a region R on the interaction graph and let E 
be the set of edges connecting sites on the boundary R C 
R to exterior sites. We are interested in the injectivity of 
the operator Tr : C^®l^l ^ C^^I^L The product form 
of the A's (HH) leads to 

{x\TR\a) = C{x) {x\F\a) , (14) 



B. Uniqueness for non-injective lattices 

Although injectivity is a generic and powerful condi- 
tion for proving uniqueness of the ground state it is not 
a necessary requirement. In fact, in the case of a PEPS 
associated to a classical model on a square lattice, al- 
though we do not have injectivity, it is still the unique 
ground state of its parent Hamiltonian . The reason 
is that we can essentially follow the main proof given 
above without using injectivity. The only thing we need 
for that is to show that the intersection properties of the 
ranges Gn in Lemma [5] remain true. This can be shown 
easily just using p4p and some dimension considerations. 
We will illustrate them in a particular case, that is, we 



6 



will prove that, with the obvious notation, 

i-i-i i-i-i i-i-i-i 

i-i-i i-i-i i-i-i-i 



(15) 



with each dot corresponding to a two-dimensional Hilbert 
space. 

Since the inclusion D is trivial it is enough to show 
that the dimension of the left hand side is < 2^° (that, 
by (|14p is the dimension of the right hand side). Now, 
the operator 1 (g) |0)Q(0|a acting on 



Hi 



Gm-m-m 

Hi 



has trivial kernel, where a is the position (2,3). To see 
this, it is enough to notice that by (fT4 



dim(l® |0)a(OU(G.-.-.)) 

1:1:1 



In particular 1 (g) |0)a(0|a has also trivial kernel when 
acting on the left hand side of (fT5| . Therefore, the di- 
mension of this left hand side is 

< dim(l ® |0)„(OU(G.-.-. ® n)) = 2^°, 

1:1:1 

where the last equality comes from the fact that, again 

by m, 

dim(l(g|0)a(0|a(G.-.-.)) = 2^. 

1:1:1 

Finally, it has to be noticed that, if we consider the 
following regions. 



Ri 



• • • 

• • • 

• • • 



i?2 



• • • 



by dimensional considerations using as above, it is 
straight forward to verify that Pr^ — Pn^ ® 1, where 
Pfl is the projection onto the orthogonal complement of 
Gr. Then, the parent Hamiltonian of our PEPS can be 
considered to be both Pr^ or Pr^ . While Pr^ has still a 
square structure, Pr^ is more natural in the sense that it 
reflects exactly the nearest neighbors of the spin in the 
center of a cross (for this reason it is the one that will 
appear in the next section). 



VII. A COMPUTABLE SUFFICIENT 
CONDITION FOR AN ENERGY GAP 

The detection of either criticality or the occurrence 
of a spectral gap is an important problem in both con- 
densed matter theory and quantum information theory. 
Unfortunately in the 2D situation there are very few cri- 
teria for the existence of a gap above the ground state 



energy. In this section we will provide a computable suf- 
ficient condition for a gap in the parent Hamiltonian of 
a PEPS. The idea comes from the ID case and the 
proof is essentially the same. For simplicity we consider 
a translational invariant frustration free local Hamilto- 
nian on a 2D square lattice. Let ft, be a locally acting 
projector and /ly its translate by i columns and j rows 
so that H = Since H is local, there exists a 

smah / C {1, . . . , iV} X {1, . . . , Af} such that hh^j > 
whenever ^ /. Then, 



Proposition 4. If H has a unique ground state (i.e., it 
is the parent Hamiltonian of an infective PEPS ) and 



hah > 



1 



1 



then there is an e > such that > eH and hence 
there is a uniform (independent of the size of the system) 
energy gap for H . 



by 



In general, one can replace the right hand side of ([T 



CHQoh 



In the case of the PEPS associated to the classi- 
cal Ising model on the square lattice, the above criterion 
finds a gap for j3 < 0.27. In this case one can indeed 
show analytically the existence of an energy gap for the 
whole range P < Pc = ^{l + V2) 0. 

To do so we need to use the concept of a Q-matrix: a 
matrix such that qa < for all i, qij > for all i j, 
and qij = for all j. As a consequence of [13, Thm 
2.1.2] the eigenvalues of Q verify = Xq > Xi > . ■ . > 
Xn-i > — 1. The Q-matrix we are going to use is the one 
that generates the Glauber dynamics associated to the 
classical Ising model That is, for a function / on the 
configuration space we define the gradient of a function 
/ at the site i as the function \Jif{x) — /(a;*) — f{x), 
where is the configuration x after flipping the spin at 
site i. With this notation we define 



where the sum is on the sites i of the lattice and the tran- 
sition rates c{i,x) are those associated to the Metropolis 
dynamics 

c{i,x) =min{l,e-''(^'-f^')(^'}, 

with Hi{x) — — XiXj (the sum on the j's connected 
by i in the interaction graph; in the case of the square 
lattice the nearest neighbors of i). 

Q is trivially a Q-matrix, and it is proven in that 
gap((5) = — Ai has a uniform (independent of the size of 
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the system) lower bound. Now we consider the Hamilto- 
nian in our quantum system given by Hq — —Q. If we 
define the matrices Qi by 

{-c{i,x), y^x 
, else 

we have that 

(i) Qi is a Q-matrix (and hence < 0), 

(ii) g = E,Q„ 

(iii) Qi acts only on the nearest neighbors of i, 

(iv) QIV') = [11, Thm. 3.5.5] and hence g.lV-) = for 
all i. 

With these properties and using it is easy to show 
that Hq is in matrix ordering upper bounded by the par- 
ent Hamiltonian of our PEPS. This immediately implies 
that the parent Hamiltonian is gapped. 

Acknowledgements: D.P.-G. was supported by the 
Spanish grants MTM2005-00082 and Ramon y Cajal. 



VIII. APPENDIX: SITE-INDEPENDENT 
TENSORS 

Here we prove that for every PEPS which is trans- 
lational invariant on a iV x M 2D square lattice there 
exists a PEPS representation with site-independent ten- 
sors. The proof follows closely the ID case 



For doing tensor contraction we can apply Dirac nota- 
tion along the edges of the lattice, but we will need some 
conventions. The first is that, when we do a horizontal 
(resp. vertical) contraction, the vertical (resp. horizon- 
tal) indices always tensorize. For example, the result 
from contracting the positions (1, 1) and (1, 2) will be 

E 4li]uudAi\rl\u2,d2\^i)h{r2\h\ui,U2)^{di,d2\v- 

li , ui , di 
r2,U2, d2 
a 

The other convention is that the edges of each row (col- 
umn) are also joined. 

We start with a general PEPS: 

l^>= E C ((A(;2 )(,,,)) |^(l,l))...|^(^,M)) 

*(l,l)i---i*(JV,M) 

and define site-independent tensors [N M)TrM Si as 
E ^uA,i,A3^ k, l)h{j,k + l, rUlfc, J, u),{k,j + 1, dU. 

u,d,l,r.j,k 

With this choice translational invariance immediately im- 
plies that 

Iv) ^ E ^ ii^H,,k))u.k)) • • • l*(iV,Af))- 

*(1,1) !■■■!*(«, M) 
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